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We present a comprehensive theory of spin transport in a non- degenerate semiconductor that is in 
contact with multiple ferromagnetic terminals. The spin dynamics in the semiconductor is studied 
during a perturbation of a general, non-collinear magnetization configuration and a method is shown 
to identify the various configurations from current signals. The conventional Landauer-Biittiker 
description for spin transport across Schottky contacts is generalized by the use of a non-linearized 
I-V relation, and it is extended by taking into account non-coherent transport mechanisms. The 
theory is used to analyze a three terminal lateral structure where a significant difference in the spin 
accumulation profile is found when comparing the results of this model with the conventional model. 

PACS numbers: 



I. INTRODUCTION 

Hybrid semiconductor /ferromagnet material systems 
play a key role in spintronics research [1]. The motiva- 
tion to study these systems is twofold. First, computing 
technologies rely on the ability to easily tune the carrier 
density in semiconductors. Second, the advances in stor- 
age applications rely on the ability to inject or extract 
spin polarized electrons across interfaces between non- 
magnetic and ferromagnetic materials [2|-0- In the last 
decade, spin injection from ferromagnetic materials into 
semiconductors has been showing a significant progress 
[7]- [17] together with a better understanding of the in- 
terface transport properties [^-[25[. This progress has 
been accompanied with theoretical analysis of basic spin 
transport phenomena starting with the conductivity mis- 
match between a magnetic metal and a semiconductor 
[2I, [26I, [27I , [28I , [2^ , and continuing with effects of electri- 
cal fields jlo, HH , of lateral transport [32[ , and of time 
dependent response [33|, [S^, l35[ . In this paper we provide 
a comprehensive theory of time-dependent spin trans- 
port in semiconductor /ferromagnet (SC/FM) systems. It 
studies the potential and spin accumulation profiles in a 
general non-collinear setup of magnetization directions. 
Two new aspects are provided in deriving the transport 
equations. First, we elaborate on the quasi-neutrality 
approximation that simplifies the description of the drift 
diffusion equations. We show that the reasoning for this 
often-invoked approximation is different than what has 
been assumed since the late 1940's [36|-[40|. Second, 
we take into consideration the localization of electrons 
due to the doping inhomogeneity of typical SC/FM junc- 
tions [1, 0, [ll| . This leads to a change of the canonical 
boundary conditions that rely on the Landauer-Biittiker 
formalism [26]- [35,]. We use this model to analyze lat- 
eral geometries which capture the vast majority of spin 
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injection experiments |4l|-[43]. The analysis considers 
the intrinsic capacitance of the SC/FM contacts and the 
two-dimensional profiles of the electrical field and spin 

accumulation [33, [si. 

This paper is organized as follows. Sec. [Ill presents a 
time dependent analysis of spin transport in a bulk semi- 
conductor region and across a SC/FM junction. Sec. Hill 
discusses the modifications that should be introduced in 
realistic systems. It deals with the revision of the bound- 
ary conditions in forward biased junctions and with their 
voltage bias limitations. In Sec. [IVlwe apply our model 
to a non-collinear, three terminal planar geometry and 
we quantify the time dependent readout across a capaci- 
tive barrier. Sec. [Vl provides a summary. Descriptions of 
technical numerical procedures are provided in separate 
appendices. 



II. GENERAL FORMALISM 

Theoretical analysis of spin injection from metals into 
semiconductors shows that non-ohmic junctions are nec- 
essary for the current to be polarized [i, [H, [53, HI, • 
More precisely, since the spin-depth conductance of the 
semiconductor (conductivity divided by spin-diffusion 
length) is much smaller than its metal counterpart, for 
spin injection to occur the junction conductance has 
to be similar or smaller than the semiconductor spin- 
depth conductance. This spin injection constraint is eas- 
ily achieved by an insulator barrier or by the naturally 
formed Schottky barrier at the interface between a metal 
and a semiconductor |50|, l5l| . In the case of a thin tun- 
neling barrier, the electrochemical potential is discontin- 
uous at the junction and as a result, the spin polariza- 
tion of the current is driven solely by the spin selective 
transmission across the junction. The much larger con- 
ductivity and spin-depth conductance in the ferromagnet 
render the spatial and spin dependence of its potential 
level negligible. Thus, in the following we describe the 
spin transport only in the bulk semiconductor region and 
across the SC/FM junction whereas the ferromagnet is 



2 



considered as a reservoir with a uniform potential level. 
We investigate in detail lateral systems that consists of 
ferromagnetic contacts on top of a non-degenerate semi- 
conductor channel. When applicable, we rely on previ- 
ous theoretical investigations of spin-transport in metals. 
These include both time dependent ^-[U and non- 
collinear [55*]- [62] aspects. In our analysis, we do not 
consider the effects of ballistic transport [59] , of external 
magnetic fields [60] , or of anisotropic spin relaxation . 



A. Bulk Semiconductor 

Macroscopic transport equations describe particle con- 
servation and current processes. These equations can be 
derived from the zeroth and first moments of the dynam- 
ical Boltzmann transport equation [H^ and they provide 
a spatial and temporal connection between spin depen- 
dent electron and current densities (n±(r,t) & j±(r, ^)). 
The accumulated (depleted) spin population at (r, t) is 
directed in the + (— ) direction. Using the relaxation 
time approximation, these derived transport equations 
are given by. 



dt 



V J, 



(1) 
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where the indices s &l s' denote either (+, — ) or (—,+). 
g > is the elementary charge and E is the macroscopic 
electric field. The spin-dependent macroscopic parame- 
ters are the spin-fiip time from spin s to spin s' and vice 
versa {ts^s' & ^s',s)^ the diffusion coefficients (Ds), the 
conductivities (erg) and the momentum relaxation times 
{Ts^m)- The current terms are eliminated by substituting 
the continuity equation into the divergence of the current 
equation. 
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At this phase, one can derive the dynamical spin depen- 
dent drift-diffusion equation by applying a series of con- 
trolled approximations after which the first line is rewrit- 
ten in a more compact form and the second line (curly 
brackets) is neglected [26|-[35|. In non-degenerate and 
homogeneous semiconductors, the diffusion constant and 
momentum relaxation time are spin and position inde- 
pendent: Dj^=D-=D & Tj^^rn='T-^m=^m' addition, 
the spin-flip times are equal and much greater than the 
momentum relaxation time, r+^_=r_^+=2rs/ r^. We 
can therefore accurately approximate the above equation 
as, 



dt 
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where u denotes the mobility (dg = qvUg). Zhu et al. 
have studied the wave-like behavior due to the second 
order time derivative in magnetic metallic systems (van- 
ishing E terms) [54]. They have shown that this effect 
becomes signiflcant at time scales shorter than r^. On 
the other hand, if the interest is in semiconductors and 
in much longer time scales then a different approach is 
needed. First, we deflne the spin polarization along the 
± axis and the charge accumulation. 
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no denotes the electron density in the conduction band 
due to the background doping. Taking into account the 
Poisson equation as well as the difference and sum of 
Eq. (|4]) with its corresponding equation (+ ^ — ) one 
gets. 
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Td is the dielectric relaxation time deflned by the ratio be- 
tween the static dielectric constant and the total conduc- 
tivity, esc/c^o, where do = qvriQ. The plasma frequency 
is deflned by ujp l/^JT^ ^J q^n^jmsc^sc where rrisc 
is the effective mass of the electron in the semiconduc- 
tor. Since in most cases the interest is in time scales much 
longer than the momentum relaxation time, it is common 
to neglect the wave-like behavior already when describ- 
ing the current components (rs^m—^ in Eq. ([2|)). The 
resulting charge dynamics is then described by a diffu- 
sion equation which in the linear regime of small charge 
perturbations reads dp/dt = —p/rd + DV^p. The ar- 
gument for invoking the quasi-neutrality approximation 
is then that any local charge imbalance (p ^ 0) is being 
screened out within a time scale of the order of Td- This is 
a widely used argument, whose origin can be traced back 
to the seminal works on bipolar transport in homoge- 
neous semiconductors [36|-0. However, by keeping the 
wave-like terms then the decay of p is actually governed 
by Tm and it is nearly independent of Td- This is a man- 
ifestation of the flnite propagation velocity which also 
results in an oscillatory behavior during the relaxation. 
It is still justifles, however, to assign p=0 in Eqs. ([6]) 
and (|8|) but the argument should refer to interest in time 
scales much longer than r^. This statement is general 
and should not change qualitatively if the original trans- 
port equations (([I]) & ([2])) are derived without employing 
the relaxation-time approximation. Moreover, the exclu- 
sion of spin dependent parameters in Eqs. ([7j) and ([8|) also 
suggests that the charge dynamics is general. For exam- 
ple, similar charge dynamics describes bipolar transport 
in homogeneous semiconductors. The charge accumula- 
tion is then due to deviation of hole and electron densities 
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from their local equilibrium, p = (An — Ap)/(no +Po)- 
In a different publication we will revisit the widely used 
quasi- neutrality concept, and elucidate the true nature of 
ultra-fast charge dynamics in various systems [64]. Here, 
we provide an example that investigates the applicability 
of the adiabatic approximations (rm = 0, p = 0) in deriv- 
ing the spin dependent transport equations in semicon- 
ductors. One should recall, however, that this classical 
approach (Eqs. ([7j) & ([8])) neglects the effect of dynam- 
ical screening. At relatively low electron densities (e.g., 
non-degenerate semiconductors) this classical description 
is accurate since the screening length is larger or compa- 
rable to the mean free path. 

We study the charge accumulation evolution p(r, t) 
after a disturbance that locally breaks the charge 
neutrality in an overall neutral bulk semiconductor. The 
charge current density at the boundaries of the system 
is Jo at all times. The dynamical response is similar 
for all systems in which the external electric field is 
negligible compared with the initial built-in electric field 
due to the charge imbalance (Jo is smaller than some 
critical current density). For simplicity the analysis 
proceeds with Jo = 0. We assume the disturbance 
happens far from the system boundaries and that it 
is spherically symmetric. Based on these conditions, 
the charge evolution and the electric field possess a 
spherical symmetry. We choose a representative initial 
charge profile due to the disturbance, p(r, t = 0) 



T .=2.75 ps 



exp(- 



C (exp(- 
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■) +exp(- 
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C is a constant chosen to keep the integrated space 
charge zero and po reflects the initial intensity of the 
disturbance at the center. The spatial extent of the 
disturbance is determined by A. For the second needed 
initial information we set dp{r,t)/dt\t=o = (this choice 
does not qualitatively change the following discussion). 

Fig. [1] shows the evolution of a charge disturbance 
whose peak intensity is po = 0.1 and its spatial ex- 
tent is A=60 nm. We consider a room temperature, 
non-degenerate n-type GaAs with a momentum relax- 
ation time of Tm = 100 fs. The resulting mobility and 
diffusion constant are, respectively, u ^i:^ 2600 cm^/V-s 
and D ^68 cm^/s. Panels (a)-(d) show, respectively, 
the evolution with these parameters for doping densi- 
ties of no - 10l^ 10l^ 10^^, & 10^^ cm-3. The resulted 
dielectric relaxation time is changed over three orders 
of magnitude {rd = Csc/ {q^tio)). In spite of the large 
changes in r^^, the decay time is about 2r^=200 fs in all 
cases. The exp(— t/2rm) decay is a universal behavior 
if 2Tm^p > 1 [64]. The results also show a clear os- 
cillatory behavior at shorter dielectric relaxation times 
where the oscillation frequency matches the plasma fre- 
quency, cOp. To understand this behavior we consider the 
Fourier transform of the initial disturbance. Its effective 
width is ~ 1 / A and its coherence time scale is defined by 
Tc = A^/D. If the initial disturbance is relatively wide 
such that, Tc » Td^ then the oscillations are governed 
by the (central) plasma frequency. The oscillatory be- 
havior is damped when Tc < Td due to the destructive 
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FIG. 1: (Color online) Charge disturbance evolution and 
propagation p{r,t). (a)-(d) are the results for different dielec- 
tric relaxation times, Td = {2.75 ps, 275 fs, 27.5 fs k 2.75 fs} 
which corresponds, respectively, to equilibrium electron den- 
sities of no = {10^^ 10^^ 10^'^ & 10^^} cm-^ The relaxation 
dynamics is similar and only slightly affected by the dielec- 
tric relaxation time. The decay time is about 2Tm=200 fs. At 
high densities, the oscillations are at the plasma frequency. 



interference between the wavevector components of the 
disturbance. Studying these and other effects (e.g., the 
role of momentum relaxation time, the non-linear terms, 
confined and open systems) are beyond the scope of this 
paper and will be studied elsewhere [64]. 

To summarize the quasi- neutrality aspect, if the in- 
terest is in spin phenomena at time scales much longer 
than Tm (and not r^^), then it is accurate to apply the 
adiabatic approximations (r^ = 0,p = 0), and get a 
divergence-free electric field and a linear dynamical spin- 
drift-diffusion equation. 
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where £sf = \/^^sj is the spin-diffusion length and the 
Einstein relation was invoked, D/v = Vt = ksT/q. 

The spin dependent electrochemical potential, is 
also an important transport quantity which has the fol- 
lowing relation with p. 



M±(r,t) Po 



-ln(l±p(r,t)) , 



(11) 



where po = Pc — <70(i*, t) denotes the spin independent 
part defined by the sum of a constant chemical potential 
and the electrical potential. The latter is driven by the 
applied bias voltage and is related to the electrical field 
via E = —Vcj). The logarithmic term refers to the non- 
degenerate case (with p = 0). 

To this point, we have treated the spin polarization 
in the channel, p, as a scalar which implicitly relies on 
the assumption that the net spin has a fixed direction 
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throughout the semiconductor channel. This description 
is vaUd in cohinear systems at which the magnetization 
directions in ah of the ferromagnetic elements share a 
common (easy) axis. In a more general, non-collinear 
configuration the boundary conditions impose a change 
in the direction of the spin polarization during the trans- 
port in the channel. For a general coordinate system 
in spin space, the spin dependent electron density is de- 
scribed by a 2x2 matrix. 



n(r,t) = ^(/ + p(r,t)-<T) , 



(12) 



where a is the Pauli matrix vector and p has the magni- 
tude p along the + direction as defined in Eq. ([5]). Using 
this notation and repeating the analysis, the components 
of the spin-drift-diffusion equation read. 




where i (j) enumerates the x, y and z coordinates in 
spin (real) space. According to Eq. ([2]), the components 
of the charge current density (vector) and of the spin 
current density (second-rank tensor) are. 



J. 
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These expressions are valid if the frequency of the applied 
electrical signal is much smaller than l/Tm- 



B. SC/FM Junction 

The description of transport is complete when the spin 
polarized currents across the SC/FM junctions are ex- 
pressed in terms of the spin polarization vector at the 
semiconductor side of the junction, p(rj,t). We follow 
the notation by Brataas et al. [59., £1], and use Eqs. (pTj) 
& (p!2|) to write the population distribution matrices on 
both sides of the junction. 



(16) 
(17) 



where e denotes the energy. As mentioned before, due 
to the conductivity mismatch the ferromagnetic side is a 
reservoir with a constant chemical potential, /io — QlV ^ 
where /io (Eq. (pT]) ) is evaluated at the semiconductor 
side of the junction and V is the voltage drop across 
the SC/FM junction. F > (V < 0) denotes forward 
(reverse) bias voltage in which electrons fiow into (from) 
the ferromagnetic contact. For compact notation, the 
boundary conditions of each SC/FM junction are writ- 
ten in a spin coordinate system at which the z-axis is 



collinear with the majority spin direction of the corre- 
sponding ferromagnetic contact. With this simplifica- 
tion, the refiection matrices are diagonal. 
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The refiection coefficients in the left (right) matrix are of 
electrons from the semiconductor (ferromagnetic) side of 
the junction. For a given material system, these coeffi- 
cients vary with the voltage drop and with the longitu- 
dinal energy, which denotes the impinging energy of 
electrons due to their motion toward the SC/FM inter- 
face. The up and down arrows denote, respectively, the 
majority and minority spin directions in the ferromag- 
netic contact where we have set the +z direction par- 
allel to the majority direction. By using the Landauer- 
Biittiker formalism, the tunnelin g cu rrent density across 
the SC/FM junction is given by [Hi, iH, El , 
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(19) 
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The first (second) term in square brackets is related to the 
transmitted current from the ferromagnet (semiconduc- 
tor) due to electrons whose total and longitudinal ener- 
gies are e and e^, respectively. The zero energy refers to 
the bottom of the semiconductor conduction band. The 
inner integration is carried over transverse wavevectors 
due to a motion in parallel to the SC/FM interface and its 
upper integration limit, /c^, denotes the wavevector am- 
plitude of an electron with energy e. In the chosen spin 
coordinate system, the transmitted spin current from the 
ferromagnetic side is non-zero only along the z direction. 
Its spin-up and spin-down components are proportional, 
respectively, to (1 — |r|p) and (1 — |r|p), where we have 



rendered the fact that 
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P. The transmit- 



ted current from the semiconductor side, on the other 
hand, includes off-diagonal mixed terms that are propor- 
tional to r^r*. This is the case when p ^ z due to the 
ffow of electrons from/into ferromagnetic contacts which 
has non-collinear magnetization directions and that are 
located within about a spin-diffusion length. By substi- 
tuting Eqs. (p!6|) - (p!8|) into Eq. ([19]) we write the energy 
resolved tunneling current density matrix. 



(1J[ 



(20) 
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The Pi components of the spin polarization vector are 
evaluated at the semiconductor side of the junction, and 
the direct and mixing conductances (per unit area) are 
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given by, 

'*-^'' = Cr^ (^^' 

In this writing, the dependence of the reflection coef- 
ficients is resolved from e & k\\. The analysis is further 
simplified if we assume that the potential level in the 
ferromagnetic contact lies beneath the edge of the semi- 
conductor conduction band, jiQ — qV < 0. Since /io < in 
a non-degenerate semiconductor , this assumption holds 
for any forward bias conditions and for relatively low re- 
verse bias conditions. Thus, we can consider only the 
Boltzmann tail of the ferromagnetic Fermi-Dirac distri- 
bution (first term in Eq. ([20]) ). In this regime, the compo- 
nents of the tunneling current density across the SC/FM 
boundaries are compactly described by, 

Ja{V) = VtG (e-^/^- - 1 - Fp,) , (23) 
Z.AV) = VtG (f (e-^/^- - l) - p,) , (24) 

^y,a{V) = 2Vt {lm[Gn\p.-MGn]Pv) > (25) 
Z.AV) = '^Vt {-MGu]P'^-MGu]Vy) ■ (26) 

The subscript a is a real space coordinate directed along 
the normal of the SC/FM interface. The total and mixing 
macroscopic conductances {G & G^) and the finesse (F) 
of the junction are given by, 

QVt Jo 
1 

Gm = -rr dee^^^-^y^-^gn{e.V), (27) 
QVt Jo 

1 

qYtC^ Jo 

The reflection coefficients that appear in the bound- 
ary currents can be either extracted from carefully de- 
signed experiments [66, 63, [1^, or calculated by various 
techniques. In this paper, we use an effective mass sin- 
gle band model [65,], [6^]-[72[. Details of the calculation 
are given in Appendix [Aj We mention that this model 
does not include the predicted effects of Ab-initio calcula- 
tions [18]- [25]. These effects may become crucial in ideal 
SC/FM interfaces and concern spin- filtering mechanisms 
or the presence of interfacial bands. 

III. REALISTIC MODELING 

There are several restrictions and mechanisms that 
should be considered in a realistic modeling of spin trans- 
port in biased SC/FM hybrid systems. The use of Pois- 
son and linear spin-drift-diffusion equations ((|9|) & (p!3|) ) 
inside the homogenous bulk semiconductor is a reliable 



macroscopic description as long as ballistic effects are 
not important. However, three crucial aspects about the 
boundary currents need to be addressed. The first re- 
lates to the inhomogeneous doping profile at the Schot- 
tky barrier and it results in a change of the boundary 
conditions. The second aspect relates to the bias voltage 
across SC/FM junctions and it limits the applicability 
of the boundary conditions to a finite bias voltage range. 
The third relates to the intrinsic capacitance of the Schot- 
tky barriers and it plays a role in the dynamics. 



A. Inhomogeneous doping profile 

At low bias voltages, the width of the barrier should be 
about or less than 10 nm in order to suppress thermionic 
currents at room temperature (see Appendix [A|). There- 
fore, the doping concentration at the Schottky junction 
should be highly degenerate, Ugb ^10^^ cm~^. If the bulk 
of the semiconductor is non-degenerate (no <dO^^ cm~^) 
then the result is a strongly inhomogeneous doping pro- 
file between these regions ^ [ol, [ll| . The need for a much 
lower carrier density in the bulk region, stems from the 
condition of optimal spin accumulation density: the bar- 
rier conductance is of the same order of magnitude as 
the semiconductor spin-depth conductance G ~ cro/hf 
[H, [32|. Since G is relatively small even with a ~10 nm 
wide barrier, the optimal 'spin-impedance' condition can 
be met with a bulk region that is moderately doped (low 
(Jo) for which the spin-diffusion length, /^j, is relatively 
large. Thus, for the spin accumulation density to be non- 
negligible the background doping densities are typically 
such that no < 10^^ cm~^. As a result of the doping in- 
homogeneity between the bulk and the barrier regions, a 
potential well is likely to be created between these regions 
^]-[75]. Even with a careful doping design at which 
there is no well in equilibrium, at forward bias when less 
electrons need to be depleted from the semiconductor, 
the well creation is inevitable. The spin related effects in 
this potential well may contribute to the spin accumula- 
tion in the bulk semiconductor region [tqI. ItII] . 

Fig. [2] shows the conduction band profile across a one- 
dimensional biased system (calculation details in Ap- 
pendix [Xj). The left and right barriers denote, respec- 
tively, the reverse (spin injection) and forward (spin ex- 
traction) biased junctions. Across the reverse biased 
junction, most of the injected hot electrons overshoot 
this region (although their reflection coefficients are slight 
modified by the potential well). On the other hand, in 
the forward direction, we should consider a new trans- 
port mechanism which involves the escape of spin polar- 
ized electrons from the potential well into the ferromag- 
netic contact. This process accompanies the previous 
mentioned process of free electrons tunneling (Eq. (p!9|)). 
Therefore, free electrons from the bulk semiconductor 
region can either tunnel directly into the ferromagnetic 
metal or feed the potential well when its localized elec- 
trons escape. To quantify the spin dependent currents 
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FIG. 2: (Color online) A scheme of the conduction band pro- 
file of a one-dimensional FM/SC/FM structure. In forward 
bias, the potential well contributes to the current across the 
interface. Inset: room-temperature J-V curve across an n- 
type GaAs/Fe junction where the semiconductor bulk doping 
is 10^^ cm""^ and the interface doping is 2 x 10^^ cm"*^. 



that flow to the potential well, we consider three mech- 
anisms. The first is the capture time of a free electron 
into the potential well. This is a spin independent time 
scale of the order of hundreds of fs and it is governed by 
spin conserving phonon-carrier or carrier-carrier scatter- 
ing processes [76|-[78|. The second mechanism is the spin 
relaxation in the potential well and its time scale is of the 
order of a few ps in III-V semiconductor quantum wells 
[79^ . The third mechanism is the spin dependent escape 
time of electrons from the i^^ localized state of the po- 
tential well into the ferromagnetic contact. Its order of 
magnitude can be calculated by a WKB method, 



1 



1 



2 Q^^sb 



■"(28) 



(29) 



05 is the difference between the work function of the 
metal and the affinity of the semiconductor (if needed it 
can also factor the pinning of the Fermi level). The bar- 
rier height from the conduction band of the semiconduc- 
tor is denoted by — Mo — qV. The localization energy 
is denoted by Ei (see Appendix [A] for its calculation) and 
uj corresponds to the parabolic curvature of the conduc- 
tion band at the barrier region. The escape time scale is 
highly sensitive to the doping level of the Schottky region. 
For example, it increases from 28 ps to 11 ns when the 
doping is reduced from Ugb ~ 2 x 10^^ to 77,55 ~ 7 x 10^^. 
These values are calculated by using GaAs bulk parame- 



ters, rrisc = 0.067mo and Cgc = 1-16 x 10~^^ F/cm and a 
typical value of (/)^ + /io — — = 0-7 eV. This process 
is spin dependent since the escape rates are proportional 
to the inverse of the wavevector in the ferromae;netic side 



(30) 



The faster escape rate of electrons with smaller wavevec- 
tor (e.g., minority electrons in iron) provides a way to 
distinguish this effect from the delocalized electron tun- 
neling whose spin polarization is opposite [i^. Due to 
the large differences between these time scales, Tcap <C 
Ts.weii <^ Tesc) wc cau assumc that (I) every electron 
that escapes from the potential well into the ferromag- 
net is being replenished immediately by an electron with 
the same spin from the bulk region and (II) inside the po- 
tential well the spin polarization is negligible, p^^ ~ 0, 
and as a result only the total current density, J^^, and 
the spin current density in the z-spin coordinate, 
are non-zero. According to our coordinate system, the 
+z coordinate denotes the majority spin direction in the 
ferromagnet. In analogy with the boundary conditions of 
delocalized electrons in Eqs. ([23 |) - ([26|) . we see that even 
when Pz=^ the spin current density in this direction is 
non-zero if F 7^0. The role of the finesse in the localized 
case is played by the spin-dependent escape times. To 
comply with this physical picture, we add phenomeno- 
logical terms to the boundary conditions of the forward 
biased junction. 



-UV) = VtG (e-^/^- - 1 - Fp,) + , (31) 
Z.AV) = FTG(F(e-^/^--l)-p,)+aI^. (32) 



The contributed current density from the potential well 
is given by. 



= 

Tp2D 



1 



^ . I esc, I 
p2Dj2D ^ 



(33) 

(34) 
(35) 



hi denotes the (bias dependent) two-dimensional density 
of electrons that are not Pauli blocked in the z*^ localized 
state. The energies of these electrons are above the po- 
tential level in the ferromagnetic contact and thus they 
can contribute to the escape process. The 1/2 factor de- 
notes the fact that the spin polarization in the potential 
well is zero (n^^^^ = ^i,-z = hi/ 2). The x and y com- 
ponents of the localized spin current densities are zero 
{piF — Py^ — 0) thus only the free electrons con- 



tribute in these directions (Eqs. ([25]) & ([26 

Incorporating the potential well contribution to the 
current density is the only way that one can fit exper- 
imental J — V curves. Using the Landauer-Biittiker for- 
malism to calculate the total tunneling current density 
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of free electrons (the trace of Eq. (p!9|) ) shows that at low 
bias conditions the current density is exponentially larger 
in the reverse direction (J(-|V|) > J{\V\) if |V|<0.2). 
This is shown at the inset of Fig.[2]for Usb = 2x10^^ cm~^ 
and no =10^^ cm~^ which are the typical experimen- 
tal values in Fe/GaAs material systems [H- lH.jis^-jiB^. 
However, the total experimentally measured J — V curves 
show that the forward bias supports larger current densi- 
ties than the reverse bias throughout this voltage range. 
Without fitting parameters, this contradiction is settled 
by including the potential well contribution [70] . There is 
an important consequence from this analysis. The escape 
current of localized electrons and the tunneling current 
of free electrons can result in opposite spin polarizations. 
Therefore, one can engineer their relative fraction by the 
doping profile [71]. For the highly doped interface, the 
escape current dominates the forward J — V curve and as 
a result the spin polarization in the semiconductor is al- 
ways along the majority spin direction (assuming that the 
semiconductor bulk region has net spin before the injec- 
tion or extraction). We will revisit this point in Sec. llVi 
In this paper, we assume that in reverse bias condi- 
tions the injected electrons overshoot the potential well. 
This approximation is accurate at very low temperatures 
where localized electrons cannot gain enough energy in 
order to pop-out of the well. However, at room temper- 
ature there is another transport mechanism to consider. 
Electrons from the ferromagnetic contact can tunnel into 
unpopulated localized states in the potential well and 
then to pop-out into the bulk region by absorption of 
a phonon or by electron-electron scattering. It should 
remain clear, however, that throughout the reverse bias 
range, the injection into free bulk states is the dominant 
mechanism due to the lower barrier that is involved in 
this process. Thus, in the following simulations we con- 
sider the localized electrons only in forward bias where 
the escape current can become the dominant transport 
mechanism. 



B. Bias voltage limitations 

Complementary magneto-optical Kerr spectroscopy 
and electrical Hanle measurements show that the spin 
polarization is appreciable only over a moderate bias volt- 
age range both in the forward and reverse directions [49| . 
In forward bias, there are two reasons that limit the 
spin polarization of extracted electrons with increasing 
the voltage across the SC/FM junction. First, at large 
positive voltages electrons can tunnel into ferromagnetic 
states above the Fermi energy where new bands with 
smaller or zero spin polarization exist [2?], [69^. This 
effect can be incorporated by calculating the reflection 
coefficients with additional bands above the Fermi level 
of the ferromagnet. Second, the barrier height is lowered 
with forward bias and thus the conductance increases 
exponentially. The spin accumulation in the semicon- 
ductor channel disappears when the barrier conductance 




FIG. 3: The ratio between the net spin after energy thermal- 
ization and the net injected spin as a function of the initial 
injected energy in a 10^^ cm~^ n-ty pe GaAs. The method of 
calculation is presented in Ref. [83] • 



significantly exceeds the spin-depth conductance of the 
semiconductor (G ^ ao/igf). At this bias regime, the 
voltage drop across the junction is negligible and as a 
result the electrochemical potential splitting between /i+ 
and /i_ becomes negligible. Both mentioned restrictions 
limit the ability to achieve spin selective extraction of 
free and localized electrons (first and second terms of 
Eq. (|32|). respectively). 

In reverse bias conditions, increasing the voltage am- 
plitude limits the spin polarization of injected electrons 
due to a transport across a wider depletion region with 
enhanced electric field [Sll, [G^, and due to an en- 
hanced spin relaxation of injected electrons prior to their 
thermalization [81, 82, 83]. The detrimental effect of the 
former can be overcome by increasing the doping con- 
centration next to the junction. However, the second 
effect is an intrinsic property that cannot be engineered 
for a given zinc-blende bulk semiconductor. Fig. [3] shows 
the fraction of the net injected spin that is left after the 
energy thermalization process as a function of the in- 
jected energy in a 10^^ cm~^ n-type GaAs [83]. The spin 
information is largely kept when the energy of injected 
electrons is less than 0.1 eV. These results are also in 
accordance with recent measured data by Crooker et al. 
[ioj . The spin relaxation of these hot electrons is gov- 
erned by the Dyakonov-Perel mechanism [sl, [85|. Due 
to the moderately reverse biased GaAs/Fe junction, we 
have assumed that the injected electrons tunnel into the 
F-valley of the conduction band. Thus, the energy ther- 
malization is governed by emission of long wavelength 
LO-phonons [8y]. This scenario is valid if the injected 
energy of hot electrons is less than 0.3 eV above the F- 
point of conduction band in GaAs. At stronger reverse 
bias conditions, the injected electrons can reach the L- 
valley and thus experience strong inter- valley scattering 
processes [s^l- In the case of a 10^^ cm~^ n-type GaAs 
at room temperature, this injection energy limit corre- 
sponds to a —0.4 V reversed biased GaAs/Fe junction 
(the Fermi energy is about 0.1 eV below the conduction 
band). 

We conclude that if the bias voltage across the SC/FM 



is moderate then one can neglect the spin relaxation pro- 
cesses during the ultra-fast thermalization of injected 
electrons to the bottom of the conduction band. This 
allows one to match the spin polarized tunneling current 
densities with the spin polarized current densities at the 
edges of the bulk semiconductor region (Eq. (p!5]) ). 

C. Intrinsic Capacitance of the Schottky barriers 

In every moment of time, the applied potential fulfills 
the Laplace equation, V^^ = V • E = with boundary 
conditions given by the charge current densities at the 
interfaces, Ja{V). In the time-dependent case the charge 
current density also includes a displacement current den- 
sity connected with charging or discharging the barrier 
capacitance (changing the width of the depletion layer), 

dV 

Jsh = CB-g^ • (36) 

V is the voltage drop across the junction and is the 
barrier capacitance per unit area whose magnitude is 
given by the ratio between the static dielectric constant 
and the width of the Schottky barrier (typically of the 
order of 10~^ F/cm^). In the following simulations, we 
will include the contribution of this current as part of the 
boundary conditions. This current can have a strong ef- 
fect on the dynamics in time scales of the order of cb/G. 

Contrary to charge currents, the spin currents are neg- 
ligibly affected by the displacement current. This is valid 
if the change in the width of the Schottky barrier, A^^, is 
such that Glsf/ (JQ ^ Ad/lsf- To derive this condition, 
we recall that the spin current densities, 3f include 
terms of the order of VrGpi. The displacement spin cur- 
rent density is proportional to qnoAddpi/dt. If we are 
interested in time scales of the order of or longer than 
the spin relaxation time {dpi/dt < pi/rsf) then one 
can readily derive the above condition. Neglecting the 
contribution of the displacement current is robust even 
at faster time scales since low voltage signals change the 
width of a heavily doped Schottky barrier by only a few 
nm (A^ <C Isf)' Barriers that are not heavily doped are 
of no interest due to the resulting negligible spin accu- 
mulation. 



IV. RESULTS AND DISCUSSION 

In showing the results of the presented theory we fo- 
cus on two aspects. First, by using a steady state anal- 
ysis we show how the spin transport in a semiconductor 
channel is affected by the escape current density and by 
the non-linear J-V relation across the SC/FM junctions. 
When contrasted with conventional spin transport anal- 
yses the results are significantly different . Second, by 
using a dynamical analysis we study the non-collinear 
magnetization configuration effect on the spin accumula- 
tion and current signals. The static and dynamical trans- 
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FIG. 4: (Color online) The lateral structure we use in the 
simulations of Sec. lIVI It consists of a semiconductor channel 
covered by three ferromagnetic terminals. The geometrical 
parameters are d—1 /xm, nm and w2—200 nm. 

In (a) the capacitor is connected to the right contact and in 
(b) to the middle contact. In both cases, the capacitance is 
C=4 fF. The magnetic configuration is denoted by the 01, 02 
and 03 angles in the ^/^-plane (measured from the +z axis). 



port analyses are performed on a lateral semiconductor 
channel covered by three ferromagnetic terminals. Fig. [H 
shows two bias settings of the studied structure where in 
each case two terminals are biased and a third terminal 
is connected in series to an external capacitor (not to be 
confused with the intrinsic Schottky capacitance of each 
SC/FM junction). By perturbing the magnetization vec- 
tor of the left or right terminals, the resulting transient 
charge current across the external capacitor allows us to 
study the spin dynamics in the semiconductor channel. 

We perform simulations for applied voltage up to 0.3 V 
with all possible in-plane magnetization alignments at 
an interval of 7r/4 in each terminal. We recall that 
the boundary conditions across a SC/FM junction were 
derived using the assumption that the spin-z axis is 
collinear with the majority spin direction in the FM (sec- 
tions [nB] and [IIlAl) • However, in order to consider all 
(non-collinear) ferromagnetic terminals and the semicon- 
ductor channel as one system, one should use a single spin 
reference coordinate system. Specifically, we transform 
the general expressions in Eqs. ([19]) & (|34|) into this new 
'contact-independent' coordinate system. The details of 
this transformation as well as the numerical procedure 
are explained in Appendix [Bl We use room tempera- 
ture GaAs/Fe material system where the semiconduc- 
tor parameters are: no=10^^ cm~^, ^^=2700 cm^/V-s, 
rs/=0.2 ns [88], msc=0.067mo. The Fermi wavevectors 
for majority and minority electrons in the iron termi- 
nals are, respectively, 1.1 A and 0.42 A where their 
mass is of free electrons [55[. The doping and static 
dielectric constant in the Schottky barrier region are 
ns5=2xl0^^ cm~^ and esc=1.16x 10~^^ F/cm, respec- 
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tively. The height of the barrier in equihbrium (V = 0) is 
0B=O-7 eV from the ferromagnet's Fermi energy. These 
barrier parameters yield a single localized energy level. 
In addition, the combined conductance (of both biased 
barriers) roughly matches the semiconductor spin-depth 
conductance {ao/isf ~4xl0^ Q~^cm~'^). This allows us 
to study cases in which the spin polarization in the chan- 
nel is relatively large {pno is of the order of no). Spin 
injection experiments, on the other hand, are currently 
limited to much smaller polarization values due to the 
self-compensation issue of silicon donors in GaAs pj| . 
This limits the effective interface donor doping levels to 
ns6~5x 10^^ cm~^ and the resulting conductance at room 
temperature to be less than 10^ rt~^cm~'^ [10]. Breaking 
this impasse is a central challenge in the way to realize 
room temperature GaAs/Fe spintronics devices. 



A. Static Results 

We use a steady state analysis to show the effects 
of the escape current density and of the non-linear J- 
V relation. First, we employ the setting of Fig. [4l^a) 
in which the right terminal is outside the path of the 
steady state charge current. The magnetization direc- 
tions of all three contacts are set parallel to the +z 
direction (^i = ^2 = ^3 = 0) so that pz{x^y) is the 
only non-zero spin polarization component. The applied 
voltage is Vdd=0.3 V for which the resulting small volt- 
age drop across each SC/FM junction justifies the use 
of our boundary conditions. Using these parameters. 
Fig. [5] shows a direct comparison of the spin polariza- 
tion in the semiconductor channel between (a) the full 
model, (b) the model without the escape current mech- 
anism and (c) the conventional model. Fig. [S^a) shows 
the result of the full model in which the boundary con- 
ditions are given by Eqs. ([I9|), ([33|) & i^. Fig. [Sfb) 
shows the case when only Eq. (p!9|) is employed. Fig.[5{c) 
shows the results of the widely used conventional model 
in which not only the escape current mechanism is ig- 
nored but a linear form of Eq. ([19]) is used. This linear 
form is the set of Eqs. ([23|) - ([26|) with the replacement of 
the e~^/^^ — 1 terms by —V/Vr- The conductance and 
finesses values in the conventional model are extracted 
around zero bias. Appendix [A] presents the calculation 
of the various reflection coefficient combinations whose 
integrations provide the conductance and mixing conduc- 
tance values. In our simulations these values are F=0.2, 
G=2Re{G|(|)}=6Im{G|(|)}=2xlO'^ Q-^cm-^. 

To understand the spin polarization in the semicon- 
ductor channel we consider two scenarios. In case (I), 
the transport across a reverse/forward biased junction is 
dominated by free/localized electrons. The net spin po- 
larization produced by a ferromagnetic contact is there- 
fore always aligned with the majority spin direction. 
Fig. [5](a) shows the results of this scenario. Case (II) 
assumes that free electrons dominate the transport re- 
gardless of the bias direction. In this case, the net spin 




X(nm) Y(nm) 

(C) 

FIG. 5: (Color online) Spin polarization in the semiconductor 
channel with the setting of Fig. SJa). The magnetization di- 
rections are all parallel to the +z axis and the applied bias is 
Vd=0.3 V. (a) is the result of the full model, (b) is the result 
without considering escape current mechanism, and (c) is the 
result of the conventional model. 



polarization produced by a reverse/forward biased junc- 
tion is aligned with the majority /minority spin direction 
of the contact. The results of such a scenario are de- 
picted in Figs. [5](b) and (c). In non-collinear configura- 
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tions, these rules can be generalized in the following way. 
In case (I) /(II), the net spin polarization vector in the 
semiconductor channel roughly points in the vector ad- 
dition/substraction of the majority spin directions of the 
biased junctions. 

We first explain the effect of the non-linear J-V rela- 
tion via a comparison between the spin polarization in 
Figs. [Ub) and (c). The most distinct difference is the 
shift in values of pz{x^y). To understand this shift we 
recall of the Fermi level positions in the semiconductor 
channel and in the ferromagnetic terminal. In a non- 
degenerate semiconductor channel the values of the pop- 
ulation distribution matrix are very small (/sc<0.05 in 
our simulated case). In very small reverse bias and in 
any forward bias the population distribution matrix of 
the FM, //m, is also very small for electrons that can 
tunnel to/from the semiconductor. The latter constraint 
is removed in higher reverse bias conditions where the 
ferromagnetic Fermi level reaches the conduction band 
edge of the semiconductor. This is the reason for the free 
electron asymmetrical current density at moderate bias 
conditions as seen from the inset of Fig. [2l The higher 
conductance of the reverse biased junction (when con- 
sidering only free electrons) dictates the sign of the spin 
polarization in Fig.^b). On the other hand, in Fig.[5](c) 
the spin polarization is nearly symmetrical about the zero 
level at the biased part of the channel. To introduce 
asymmetry in the conventional model one has to artifi- 
cially plug different conductance values for the forward 
and reverse biased junctions. However, due to the linear 
approximation around zero-bias, the formal derivation of 
the boundary conditions in the conventional model re- 
sults in identical conductance values. 

Fig. EJa) highlights the unique behavior of the poten- 
tial well in the interface doping area. We find the oppo- 
site of the conventional result which states that antipar- 
allel configurations lead to a much larpr spin polariza- 
tion than parallel configurations [28, 89]. In collinear 
two-terminal systems the potential well flips the role 
of the parallel and antiparallel magnetization configura- 
tions (if we ignore the electrical field effect and assume 
that the channel length is smaller than the spin diffu- 
sion length.) The potential well significantly changes the 
shape of the spin accumulation where spins diffuse in 
opposite directions. This is seen by the opposite slopes 
of the spin polarization in the left part of the channel 
in Figs. [5](a) and (b). Notably, Pz has a high plateau 
below the middle forward biased contact. Since the con- 
ductance of the forward biased junction increases when 
we incorporate the escape current mechanism, the rela- 
tive portion of the voltage drop across the semiconductor 
channel increases and as a result the electric field in the 
channel increases. The field pushes spins carried by in- 
jected electrons to the — E direction and thus spins accu- 
mulate near the forward biased side. Studying this high 
spin-polarization regime is possible when working with 
p rather than with its direction and with /i-t (Eq. (pT]) ). 
Fig. ini^a) shows the x-averaged electric potential along 



the semiconductor channel, l/hj^ dx(j){x^y). From the 
slope of the curves we can estimate the electrical field in 
the left part of the channel to be of the order of 2 kV/cm 
when incorporating the escape current mechanism and 
about four time smaller in the other cases. The resulting 
drift velocity (i^|£^|) is still below the saturation velocity. 

To explain the effect of the electric field on the spin 
currents we employ the full model with the setting of 
Fig.m^b) and Vdd=0.l V. Note that both charge and spin 
currents can flow in the semiconductor channel under the 
floating middle contact. Figs. [6l[b)-(d) show the aver- 
aged spin components along the semiconductor channel 
{1/h Jq^ dxpi{x, y)). As can be seen from Eq. ([15]), the sig- 
nature of the electric field is evident at fields amplitudes 
which exceed Vr/isf |3Q| (>250 V/cm in non-degenerate 
n-type GaAs at room temperature). Fig. [6l^b) shows the 
averaged Pz component for five magnetization directions 
of the right biased contact whereas the other two contacts 
are set parallel in the +z direction (^i =03 =0). The 
electric field opposes the diffusion of spins away from the 
forward biased junction and the spin polarization from 
the reverse biased junction spreads throughout the chan- 
nel. This is the reason that even in the antiparallel con- 
figuration (02 = tt), Pz is not much into negative values 
beneath the forward biased (right) contact. 

Figs, injc) and (d) show, respectively, the averaged Py 
and Px for 0i =02 =0. Here the non-collinearity of 
the floating middle contact 03 (from 7r/2 to Sir/ 2) is the 
only drive for the Py and Px components. As before, 
the spin accumulation is pushed by the strong electric 
field toward the forward biased junction in the right part 
of the channel. The out-of-plane spin component is a 
useful probe of spin polarization beneath ferromagnetic 
contacts. For the chosen setting, we note that Py which 
is a direct result of the non-collinear middle contact is 
smaller than px. The reason is that Px is a mixed term 
that is proportional to cross product of the spin polariza- 
tion vector and the magnetization direction in the middle 
contact (p X 1113). The spin-polarization in the channel 
is nearly collinear with the z-axis because of the magne- 
tization directions of the left and right biased contacts 
(01=03=0). Taking, as an example, the case that the 
magnetization in the floating middle contact is along the 
?/-axis (03=7r/2), then the mixed Px component can be 
larger than py due to the relatively large Pz in the chan- 
nel. On the other hand, the Py component is generated 
by the very small voltage, V/, across the floating con- 
tact. This small bias results in a charge current density 
of magnitude GVj which due to the external capacitor 
is contrasted by an equivalent charge current density of 
magnitude GVtF{p • 1113) (nullifying the expression on 
the right hand side of Eq. ([23]) with V = Vf ^0). 

B. Dynamic Results 

We use a dynamical analysis to study how the non- 
collinearity affects the current signals. The setting of 
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FIG. 6: (Color online) (a) The averaged electrochemical po- 
tential along the semiconductor channel for the setting of 
Fig. HJa) where ^^^=0.3 V and 0i=02=03=O. (b)-(d) The 
components of the averaged spin polarization vector along 
the semiconductor channel for various magnetization config- 
urations with the setting of Fig. IHb) where Vdd=0.1 V and 
01=0. In (b) the middle floating contact is also fixed at 03=0. 
In (c) and (d) the left biased contact is fixed at 02=0. 



Fig.m^a) is used where the magnetization direction of the 
middle terminal is fixed along +z (02=0). The right ter- 
minal is outside the path of the steady state charge cur- 
rent, and its magnetization direction is perturbed accord- 
ing to 03 (t) = cos(27rt/rr)i + sin(27rt/rr)^ for < t < Tr. 
The transient current across the external capacitor is 
then evaluated for various magnetization directions of the 
left terminal (0i). Similar dynamical setups have been 
suggested for collinear configuration by Cywinski et al. 
[34]. Here, we allow for a more flexible operation regime 
(non-collinearity) and we offer a different physical inter- 
pretation of the dynamical results while emphasizing the 
robustness of the signals' signature. Fig. [71 shows the 
transient currents across the external capacitor for four 
01 cases. The applied bias is Vd =0.1 V, the external 
capacitance is C = 4 £F, and the depth of the system 
in z direction is 1 /im. The magnetization direction of 
the right terminal completes a single clockwise rotation 
in rr=3 ns. 

The transient currents in Fig. [71 are described by 
— {C/q)diir/dt where /ir(^) is the potential level in the 
right terminal. This current is also the integrated cur- 
rent density at the right contact. This current density is 
given by. 



JrH) = - \g 



n2D 



(1 



-GVTFp,{t)+Cb 



dV{t) 



dt 



e(FW)j vit) 

(37) 



the right terminal, and Tesc,o is the escape time at volt- 
age drop. The terms that involve G are from lineariz- 
ing Eq. ([23]) around V=0 and the C5 term is due to the 
intrinsic capacitance of the Schottky barrier (Eq. ([36]); 
we use C5=10~^ F/cm^ in the simulations). The term 
that involves the step function, 0(V(t)), is due to the 
escape of localized electrons (linearizing Eq. ([33]) around 
V=0). riw denotes the electron density in the poten- 
tial well and n2D = qVTrrisc/^^^ • As discussed at the 
end of Sec. IIII A[ our modeling includes transport of lo- 
calized electrons only at positive voltages. As a result, 
the effective barrier conductance is discontinuous at zero 
bias. This discontinuity is the reason for the 'cusp' points 
at times smaller than 3 ns in Fig. [3 Since we neglect 
the transport mechanisms that involve the potential well 
when V <0, this discontinuity is a model dependent ar- 
tifact. Pz{t) is the projected spin polarization vector on 
03 (t). The spin polarization vector is nearly constant 
and points in the vector addition of the majority spin di- 
rections of the biased ferromagnetic terminals (without 
considering the escape currents it would be the vector 
substraction). Thus, 



cos(^)cos(^ 



<t <rr 



3 ns 



otherwise 



(38) 



V{t) denotes the (small) self- adjusted voltage drop across 



For cases that 0i 7^ {0, tt}, the discontinuity in dpz{t)/dt 
at t=rr results in an additional 'cusp' point at this time 
(see Fig. [7] at t=3 ns). At times greater than 3 ns, Jr 
is governed solely by the dynamics of V{t) towards its 
original value prior to the perturbation. At shorter times, 
Jr is governed by the counteracting response of V{t) to 
the perturbing Pz{t). This response aims at finding a 
new steady-state condition and its delay time is dictated 
by simple circuit analysis (see Fig. 1 in Ref. [13] )• If the 
delay is much longer than the rotation time then V{t) can 
be viewed as static. In this case, by inspection of Eq. ([37j) 
we see that JR{t) follows the shape of Pz{t) and its peak 
reaches an optimal value of VtGF\p\ (independent of the 
capacitance). The drawback of using a long delay is due 
to the slow dynamics at t> Tr. On the other hand, if the 
delay is very short then V{t) adiabatically follows Pz{t)- 
However, the resulting peak is now smaller (roughly) by 
the ratio of the delay and rotation times. 

We use the above analysis to elucidate some of the 
general features of the current signals by concentrating 
on the 01 = case of Fig. [71 This curve shows two 
'cusp' points around 1 ns and 2.35 ns. These are the 
times at which V{t) changes its sign. If the response 
time of the system was instantaneous (zero capacitance) 
then V{t) would have followed (^) without a delay and 
these points would appear in 0.75 ns and 2.25 ns (where 
Pz{t) changes its sign). The reason for the longer de- 
lay in the first point (1 ns - 0.75 ns) than in the second 
point (2.35 ns - 2.25 ns) is due to the larger effective 
barrier conductance in forward bias. The initial current 
shape between to 1 ns is due to the initial decrease 
of Pz{t) (Eq. (|38l) with 0i=O). This change is counter- 
acted by a 0.25 ns delayed increase of V{t) that tries 
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FIG. 7: (Color online) Transient currents across the external 
capacitor due to a single 3 ns clockwise rotation of the magne- 
tization direction of the right terminal. The inset shows the 
times at which current changes sign from positive to negative 
as a function of 01. 



FIG. 8: (Color online) Transient currents across the external 
capacitor due to a single 3 ns clockwise rotation of the mag- 
netization direction of the left terminal. The inset shows the 
times at which current changes sign from negative to positive 
as a function of 03. 



to establish a new steady state. In the second branch, 
1 ns< t <2.35 ns, V{t) is positive in order to counteract 
the negative Pz{t). The total transient current begin to 
decrease due to the turn-on of the escape current process. 
In the third branch, 2.35 ns< t <3 ns, Pz{t) is positive 
again and the sudden current drop at 2.35 ns is due to 
the stop of the escape current (the —GVrFpzit) com- 
ponent of the current is counteracted by a weaker and 
slower response of V{t)). One can repeat this analysis 
for the other signals in Fig. [71 The difference in their 
shapes is governed by the 0i phase term of Pz{t). As a 
result, these signals have an apparent trend in shifting 
the time at which the current changes sign from positive 
to negative. This crossing time, denoted by to, shows a 
linear dependence in 0i as can be seen from the inset of 
Fig.E 

Fig. [8] shows the current signals for the opposite case 
in which the left magnetization is rotated and the right 
magnetization is set at various directions (with the same 
bias setting as before). We observe similar patterns. In 
this operational regime, the spin polarization vector p 
changes in time beneath the right terminal whereas 03 
is constant. The response time is longer since the infor- 
mation needs to pass the delay of two (rather than one) 
Schottky barriers. For this setting, we denote to as the 
time at which the signal switches sign from negative to 
positive. We observe a similar and nearly linear rela- 
tion between 03 and to. However, the slope of the line is 
about twice as much then before. This double spacing is 
best explained with the vector addition rule we have pro- 
vided in the static regime. The spin polarization vector 
in the semiconductor channel points roughly to the mid- 
way between 02 = ^ and 0i. Thus, the rotation speed 
of p beneath the floating contact is about half that of 
01. In terms of distinguishing different states, the sec- 
ond setting has a doubled time resolution compared with 
the previous case. If we employ ^^^=0-3 V, then the pat- 
terns of the signal in each setting are very close to the 



above cases, while the scale of the signal is about 3 to 
5 times larger and the delay between Pz{t) and V{t) is 
shorter. Because of its physical origin, the current signal 
patterns are very robust and universal. Discussion about 
the coercivity and noise concerns can be found in Ref. 

C. Configuration Analysis 

In this part we explore the maximal number of config- 
urations that can possibly be stored in a system (as well 
as how to implement them) . When two biased ferromag- 
netic contacts have a certain magnetization alignment, 
the spin polarization vector results in a unique setting 
in the semiconductor channel. The vector is more or 
less constant throughout the channel if the spin diffusion 
length is longer than the channel. Thus, this vector labels 
the particular magnetization configuration. One way to 
gain access to this vector information is to measure the 
total resistance of the two terminal system [59]. This 
measurement can tell the relative angle but not the mag- 
netization direction of each contact. In order to extract 
the information stored in these magnetization directions 
completely, we add a third ferromagnetic contact, float- 
ing or semi- floating, on top of the same semiconductor 
channel. We use the setting of Fig. HJa) where the third 
magnetization direction provides a reference direction. 

We have shown the physical connection between the 
measured signal and its magnetization configuration in 
the previous part, but we still need a systematical anal- 
ysis to quantify the high density storage capacity in the 
designated circuit. This analysis is based on symmetry 
considerations. We will assume that in each contact there 
are two mutually perpendicular easy axes, determined by 
the semiconductor crystal orientation and by the thick- 
ness and shape of the ferromagnetic contact [90]. Fig. 
[9] shows all of the possible magnetization configurations. 
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FIG. 9: (Color online) All possible magnetization configura- 
tions with 03 = as a reference direction. The easy (in plane) 
axes of the system are in the y and z directions. There are 
10 independent configurations in the static case. Each of the 
'off-diagonal' configurations has a respective symmetric con- 
figuration that provides identical voltage drop across the right 
contact. 

The right magnetization is unchanged due to the system 
rotational invariance to voltage or current measurement 
(03 serves as a reference direction). In principle, there 
are 10 distinct voltage output values for these 16 con- 
figurations (the output is the voltage difference between 
the floating contact and the ground). It is less than 16, 
because the voltage of the floating contact is decided by 
Pz and /io beneath it and thus the voltage has another in- 
variance when the magnetization alignment is symmetric 
about z. The four configurations on the 'diagonal' of the 
4x4 table in Fig. [9] have a unique output, whereas the 
'off-diagonal' configurations are symmetric. The degen- 
eracy of the voltage signal of all six 'off-diagonal' pairs 
can be lifted by fixing 0i and 02 while turning the 03 
direction by 7r/2 and repeating the measurement. Each 
configuration then has a unique combination of two re- 
sults. In fact, performing two static measurements for 
different 03 directions is equivalent to rotating 03 in a 
given direction and measuring the transient current dy- 
namically, as showed in the first setting in Sec. lIVBi 
In the dynamical method, we can completely distinguish 
all 16 states with a single measurement. The given ro- 
tation direction serves to break the symmetry about z. 
The dynamic determination scheme is suitable for high 
frequency operation regime by properly selecting the ex- 
ternal capacitor such that the signals are less affected by 
the noise [H,^. 

V. SUMMARY AND OUTLOOK 

We have presented a detailed model that describes the 
spin transport in hybrid semiconductor /ferromagnet sys- 
tems. The derivation of the transport equations and their 
boundary conditions are then used to model lateral hy- 
brid systems with general non-collinear magnetization 
configurations. We have corrected the arguments that 
lead to the use of the quasi-neutrality approximation and 
explain their momentum and dielectric relaxation times 
dependence. The spin currents due to tunneling of free 
electrons across a semiconductor /ferromagnet junction 
were derived using a rigorous non-linear bias dependence 



of the tunneling current. The bias voltage limitations of 
spin injection and extraction were also discussed. 

We have introduced an important contribution to the 
junction tunneling current that is governed by the es- 
cape of localized electrons. This process is a result of the 
usually employed inhomogeneous doping at the vicinity 
of the semiconductor /ferromagnet interface. The escape 
current is incorporated to the boundary conditions in for- 
ward bias. If the doping inhomogeneity is large then 
this process becomes the dominant current mechanism 
and as a result the spin accumulation patterns of parallel 
and anti-parallel configurations are flipped (compared to 
the case that only delocalized electrons are considered). 
We have provided simple rules for estimating the direc- 
tion and magnitude of the spin-polarization vector in a 
semiconductor channel that is covered by spin selective 
and biased ferromagnetic contacts. Our results illustrate 
the importance of using impedance matched tunneling 
barriers with the semiconductor spin-depth conductance. 
This matching condition enables a large spin polarization 
even if the spin-selective barriers are not ideal. We have 
introduced a dynamical method that can clearly iden- 
tify the non-collinear magnetization configuration from a 
three-terminal lateral structure. The amplitude and pat- 
tern of the current signals were explained using a spin 
dependent circuit analysis that incorporates the capaci- 
tive nature of the semiconductor /ferromagnet junction. 
The presented dynamical method can be used in spin- 
tronics devices for storage beyond the binary limit. 

This work is supported by DOD/AF/AFOSR FA9550- 
09-1-0493 and by NSF ECCS-0824075. 



APPENDIX A: BARRIER STRUCTURE DETAILS 

This appendix elaborates on the model that we use 
to calculate the localization energy {Ei) and the energy 
resolved reflection coefficients (r^ & r^). These parame- 
ters depend on the voltage and are needed for calculat- 
ing the spin dependent direct and mixing conductances 
(Eqs. ([2T]) & ([22]) ), as well as the escape times from the 
potential well (Eq. ([28]) ). 

Fig. [TOlf a) shows the energy profile that we use to cal- 
culate the localization energies and the spin-dependent 
reflection coefficients. V is the voltage drop across the 
junction. The conduction band energy in the bulk semi- 
conductor is the reference level (£^c=0). The Schottky 
barrier height is 0^=0.7 eV from the Fe Fermi energy 
which is £1=4.5 eV (6:|=0.67 eV) for majority (minor- 
ity) electrons from the bottom of the ferromagnetic con- 
duction band [55]. In the semiconductor, the doping 
in the bulk and barrier regions are no =10^^ cm~^ and 
ns5=2xl0^^ cm~^, respectively. The effective masses are 
0.067mo in GaAs and mo in Fe. The Schottky barrier is 
located in the (0, ^i) region where the conduction band is 
parabolic (obeying the Poisson equation). The doping in- 
homogeneity between the bulk and the interface regions 
generates a potential well next to the barrier [73|]-l75[. 
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FIG. 10: (Color online) (a) The energy profile of the GaAs/Fe 
junction (see text for detailed parameters), (b)-(d), 1 — |r|(|)|^ 
and Im[l — r|r|] versus the longitudinal kinetic energy for five 
voltage levels. 



We model this well by a flat potential region from li to 
^2 and then by a gradual linear increase to the bulk level 
from ^2 to ^3. The width of the potential well in its flat 
region, ^2—^1, is governed by the voltage drop across the 
junction (V). It shrinks/expands with increasing \V\ in 
reverse/forward bias conditions. This flat region vanishes 
at a reverse bias of -0.2 V where £2 = ^1 =8.5 nm. The 
depth of the potential well with respect to the conduction 
band of the bulk region is /c^Tln(ns5/no). The gradual 
region width is -^3 — ^2=4 nm. The geometrical details of 
this approximated profile do not substantially change the 
results of a rigorous self-consistent Schrodinger-Poisson 
equation set [lO, [lli • 

The spin dependence of the reflection coefficients is 
governed by the different Fermi velocities in the ferro- 
magnetic region. We calculate these coefficients by as- 
suming specular transport across the junction and apply- 
ing a transfer-matrix method for the above energy pro- 
file (bold line in Fig. HOja)). Figs. [iQKb)-(d) show the 
spin dependent transmission coefficients and the imagi- 
nary part of the mixing term as a function of the longi- 
tudinal energy (kinetic energy along the GaAs/Fe inter- 
face normal). The real part of the mixing term satisfies 
2Re[l — r^r*] ~ (2 — \r^\^ — |r| , in agreement with the 
conclusion in Ref. ^] . The potential well leads to a rel- 
atively strong transmission of low energy free electrons 
(Ramsauer-Townsend resonance). This is shown by the 
peak at low ex regime in each of the F > curves. The 
summation over the kinetic energy in the parallel plane 
smears this peak in the I-V curve (inset of Fig. [2]). 

To calculate the localization energy we solve the 



Schrodinger equation for the semiconductor part of the 
energy profile (and replace the ferromagnetic part with 
an infinitely thick barrier whose height is ^b). This pro- 
file yields a single localized energy level in V g[-0.2 V, 
0.2 V] and its value is essentially linear from -0.096 eV to 
-0.13 eV. We then calculate Tesc (Eq. ([28|) ) and hi from 

Finally, we present an analytical model for the case of 
a simple rectangular barrier. This model is then used to 
extract the direct and mixing conductance values. The 
reflection coefficient of a rectangular barrier with width 

and with barrier height (j)B from the ferromagnetic 
potential level (/io — qV) is given by [65], 



e'"^\*(^)7.c-7T(i)7.*c 



(Al) 



2ms 



i^b + iksc , 7T(i) = f^b^i 



The wavevectors ksc & ^fm,m) ^re along the normal of 
the SC/FM interface, rrisc and rrifm are, respectively, 
the electron effective mass in the semiconductor and fer- 
romagnet. The spin selectivity of the reflection is solely 
due to the spin dependent wavevectors in the ferromag- 
netic contact, %m,T(i)- 

Spin injection is possible when two conditions are met. 



L 
II. 



e''°^ > 1 , 



2mscd , _ 
A = —^kBT< 1. 



(A2) 
(A3) 



where hzo = hZb{£± =0)- In triangular or parabolic shapes 
we render the same conditions since the numerical values 
of effective tvo are somewhat less but of the same order. 
The first condition guarantees a resistive rather than an 
ohmic contact [26', l27l |. The second condition guaran- 
tees that tunneling is the dominant transport mechanism 
across the barrier whereas the thermionic current is neg- 
ligible. By increasing the electrons energy, condition II 
makes the Boltzmann tail of the population distribution 
decay faster than the increase in the transmission coef- 
ficient, 1 — Thus, the main contribution to the 
current is from electrons whose energy is at the bottom 
of the conduction band. These conditions may also be 
used to simplify the calculation of the macroscopic con- 
ductances and finesse in rectangular barriers. 



G 



F = 



Vi 



hAo yi + 1 + 7?jy ' 



(A4) 

(A5) 
(A6) 



where one can see that 2Re[G|J = G, and ^0 is defined 
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by, 

Ao (7rcZ)3/2 A-3/2 - 1 ' 

L.e=0 

^T(i) = — T— • 

fflfm riO 

In semiconductors with a small electron effective mass 
and in low bias voltages the finesse is approximately 

APPENDIX B: NUMERICAL ANALYSIS 

In this appendix, we discuss the techniques we use to 
solve the overall set of equations and we provide details 
about the numerical scheme. Essentially, we need to solve 
a system of four equations (Eq. ([9|) & Eq. ([13])) with 
proper boundary and initial conditions. Boundary con- 
ditions for both the static and dynamic cases consist of 
the continuity of charge and spin current density. We 
render the no n- approximated form of the free electrons 
boundary current densities (Eqs. (p!9|) - ([22|) ). 

We discuss the static case first. When the system is 
in a steady state, the charge and spin current densities 
(Eqs. ([Hj) & (p!5|)) have no normal components at the 
boundaries that are not in contact with ferromagnetic 
terminals. In addition, the external capacitor forces a 
zero total charge current across the terminal that is at- 
tached to it. For boundaries with biased ferromagnetic 
terminals, the normal current density component is equal 
to the corresponding current density through the SC/FM 
junction. These current densities have been derived in 
Sec. IIIBI and Sec. IIII Al under the assumption that the 
z-axis is collinear with the respective majority spin di- 
rection. However, in order to incorporate the three non- 
collinear ferromagnetic terminals and the semiconduc- 
tor channel in a single system, we need to work with a 
contact-independent reference coordinate system. Specif- 
ically, we need to transform the expression in Eqs. ([19]) 
and ([34]) into the new coordinate system. We rewrite x, 
y and z in these equations as x , ^ and z to represent 
the contact-dependent coordinate. We reserve x, y and z 
for coordinates in the contact-independent system. The 



current density expressions in these two sets of frames 
are related by 

Zy,a = M COS^ Sin^ Zy,a , (Bl) 

Zz,a ) \0 -Sin^ COS^ ) \Zz,a ) 

Note that the charge current density does not depend on 
the spin space coordinate. The components of the spin 
polarization also need to be expressed in terms of the 
components in the contact-independent system (by the 
reverse rotation transformation), 

Px\ (10 \ / 

Py = COS0 — sin^ py . (B2) 
Ps J y sin (/) coscj) ) \pz ) 
For time dependent simulations, we adopt a similar 
process except that we need to consider the displacement 
current due to capacitance embedded in the system. This 
includes the intrinsic capacitance of the Schottky barrier 
(Sec. IIII Cp and the external capacitor between ground 
and the semi-floating right terminal. The initial condi- 
tion in the dynamical case is the steady state spin polar- 
ization and electrochemical potentials of its correspond- 
ing initial configuration. 

We employ a finite difference method to obtain the spin 
polarization vector and the electrochemical potential in 
our multi-terminal system. The computational grid rep- 
resenting the two-dimensional semiconductor region has 
21 X 141 nodes with a 5 nm interval. A system of 4 
differential equations is to be solved in this region. To 
obtain the desired accuracy, we use two major steps with 
iteration methods [9l|, |92[. Step I includes the evalua- 
tion of nonlinear coefficients using the under-relaxation 
method. These coefficients are updated periodically. 
Step II solves the linearized system of equations using an 
iterative technique, where we adopt the successive-over- 
relaxation method. For time dependent simulations, we 
generalize this procedure by using a Crank-Nicolson im- 
plicit method with a time interval of 0.02 ns. The space 
and time intervals are adjustable over several orders of 
magnitude, and we choose them in consideration of both 
the result details and the computation time. 
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